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ABSTRACT 

We propose a method that can determine the ellipticities of source stars of microlensing events produced by 
binary lenses. The method is based on the fact that the products of the caustic-crossing timescale, Af, and 
the cosine of the caustic incidence angle of the source trajectory, k, of the individual caustic crossings are 
different for events involving an elliptical source, while the products are the same for events associated with a 
circular source. The product Afj^ = Af cosk corresponds to the caustic-crossing timescale when the incidence 
angle of the source trajectory is k = 0. For the unique determination of the source ellipticity, resolutions 
of at least three caustic crossings are required. Although this requirement is difficult to achieve under the 
current observational setup based on alert/follow-up mode, it will be possible with the advent of future lensing 
experiments that will survey wide fields continuously at high cadence. For typical Galactic bulge events, the 
difference in At± between caustic crossings is of the order of minutes depending on the source orientations 
and ellipticities. Considering the monitoring frequency of the future lensing surveys of ^ 6 times/hr and the 
improved photometry especially of the proposed space-based survey, we predict that ellipticity determinations 
by the proposed method will be possible for a significant fraction of multiple caustic-crossing binary lens events 
involving source stars having non-negligible ellipticities. 
Subject headings: gravitational lensing - stars: fundamental parameters 



1. EVITRODUCTION 

Microlensing was originally proposed to search for Galac- 
tic dark matter in the form of massive compact halo objects 
(|Paczvhski 1986). However, with the ability to detect a large 
number of events, especially toward the Galactic bulge field, 
and to intensively monitor the detected events from follow-up 
observations, microlensing has been developed into a power- 
ful tool to study various aspects of stellar astrophysics JGouldl 
I2OOI). Especially, intensive observations of caustic-crossing 
binary lens events make it possible to measure the limb- 
darkening structure of stars, which in essence maps the at- 
mospheric t emperature as a fu r iction of depth. 

Recently, iRattenburv et alJ J2005h demonstrated an addi- 
tional applicability of microlensing in stellar astrophysics by 
measuring the shape of the source star in microlensing event 
MOA 2002-BLG-33. Due to rotation, stars are not exactly 
circular. For example, Altair (a Aquillae) has an ellipticity 
of e = (l-b^/a^y^^ = 0.23 ("van Belle'2001'i and the Achernar 
(a Eridani), the flattest s tar ever measured, has an ellipticity 
of e = 0.59 (Domiciano de Souza et al. 2003). The star mea- 
sured bv iRattenburv et al. (2005) has an ellipticity of e = 0.2. 
Here a and b represent the angular sizes of the semi-major 
and semi-minor axes, respectively. The nature of the source 
flattening makes the lensing lightcurve differ from that of a 
circular case. However, detecting the deviation induced by 
the source flattening is very difficult because the deviation is 
much smaller than the deviation caused by limb darkening 
(see § 3). For the case of MOA 2002-BLG-33, an eflipticity 
measurement was possible because of the very special lens- 
source geometry in which the source star was closely bounded 
on all sides by the caustics of the lens. However, such a case 
is extremely rare. 



In this paper, we propose a more general method that can 
determine the source ellipticity. The method requires to re- 
solve at least three caustic crossings and precisely measure 
the durations of the individual caustic crossings. Although 
this requirement is difficult to achieve under the current ob- 
servational setup based on alert/follow-up, it can be achieved 
by future lensing experiments that can survey wide fields con- 
tinuously at high cadence. 

The paper is organized as follows. In § 2, we describe 
the basics of caustic-crossing binary lens events involving 
elliptical source stars. In § 3, we demonstrate the difficul- 
ties in noticing the elliptical nature of source stars from the 
lightcurve of a single caustic crossing. In § 4, we explain how 
the source ellipticity is measured by the proposed method and 
demonstrate the feasibility of the method. In § 5, we summa- 
rize the results and conclude. 

2. CAUSTIC crossing OF AN ELLIPTICAL SOURCE 

If a microlensing event is caused by a lens system com- 
posed of two masses, the resulting lightcurve can differ dra- 
matically from the smooth and symmetric curve of a single 
lens event. The main new feature of binary lensing is the caus- 
tics, which represent the positions in the source plane where 
the lensing magnification of a point source becomes infinite. 
The set of caustics forms one or multiple closed figures where 
each figure is composed of concave line segments (fold caus- 
tics) that meet at points (cusps). The number of separate sets 
of caustics is unity when the separation between the lens com- 
ponents is comparable to the Einstein ring radius 9e, and it 
becomes two and three when the separation is substantially 
larger and smaller than 9e, respectively. Due to the larger 
cross-section of the fold caustic compared to the cusp, major- 
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ity of caustic crossing events are fold-caustic crossings. Near 
a fold caustic, the magnification of the source star brightness 

-1 /2 

is A oc M^' , where u± is the angular normal distance of the 
source from the fold caustic in units of the Einstein ring ra- 
dius. Because of the divergent nature of the lensing magnifi- 
cation near the caustic, i.e. A — > oo as mj^ 0, the lightcurve 
of a caustic-crossing event is characterized by sharp spikes. 
Since the caustics form closed figures, caustic crossings oc- 
cur in multiple times of two (at the entrance and exit of the 
caustic). Mao & Paczvriski ( 1991) predicted that ^ 1% of all 
Galatic microlensing events would be caustic-crossing events 
and this ra te is roughly consistent with t hat detected by the 
MACHO ( Alcocketal.ll200(l) and OGLE jJaroszvriskl200l 
lensing surveys. 

The usefulness of caustic-crossing events lies in the diver- 
gent nature and large gradient of the lensing magnification 
near the caustic. When a source approaches the caustic, differ- 
ent parts of the source star surface are magnified by different 
amounts (finite-source effect). Then, the lensing magnifica- 
tion of a finite source is equivalent to the intensity-weighted 
mean value averaged over the source star surface, i.e. 



Afs(0 = 



(1) 



where A denotes the point-source magnification, C is the vec- 
tor position of the center of the source, C' is the displace- 
ment vector of a point on the source star surface with re- 
spect to the source star's center, /(C') is the intensity profile 
of the source star, and the two-dimensional integral is over 
the source-star surface S. Precise photometry during the caus- 
tic crossing enables one to mea sure the limb-darkening pro- 
files of source stars ^Witt"199j) and the limb-darkening co- 
efficients of six stars with vario us stellar types were actually 
measure by this method ( Albrow et al. 1999a b, 2000, 200 ll 

lAfonso et al 2000: Fields et al. 2003: Kubas et al. 2005). It 

is also possible to stud y irregular s tructures on the source star 
surface, such as spots jHan. et alJl2000t iChang & Hani.2002t 
[Hendry et al. 20^. 

The lensing magnification of a caustic-crossing binary lens 
event involving a circular source is characterized by eight pa- 
rameters. Three of these are identical to those of a single 
point-lens event: the Einstein timescale, Ie, the time of the 
closest source approach to the center of the binary lens sys- 
tem, to, and the separation between the source and center at 
that time, mq. For a binary lens, the lensing magnification is 
not radially symmetric and thus an additional parameter of the 
angle between the source trajectory and the binary axis, a, is 
required to characterize the source trajectory. Another two 
are the parameters related to the binary lens system: the mass 
ratio, q, and the separation (normalized by 9e), s, between 
the lens components. The final two parameters are related to 
the source star: the normalized source radius, = Oj^/Oe, and 
the limb-darkening coefficient. Here 0* represents the angu- 
lar radius of the source star For a limb-darkened source, the 
surface brightness profile is often modeled as 



m 

7(0) 



= i-r> 




(2) 



where Tx is the limb-darkening parameter and 9 is the angular 
distance measured on the source star surface from the center 
to the position within the source. In the standard formalism. 




Fig. 1. — Definitions of the source orientation angle ip, caustic incidence 
angle k, normaUzed semi-major axis pa, and the projected width pj_ of the 
source star. 



the linear limb darkening is expressed as 

1- 
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= 1 - ca 

1(0) 



1-^ 



(3) 



Then, the relationship between the expressions of the linear 
limb-darkening coefficients is 



c\ = 



3Tx 
Tx + 2' 



(4) 



For an event involving an elliptical source, two additional pa- 
rameters are needed. These are the eccentricity of the source, 
e, and the orientation angle of the source major axis with re- 
spect to the perpendicular to the caustic, t/j (see Fig.fQ. There- 
fore, the total number of parameters required for the full char- 
acterization of a caustic-crossing binary lens event associated 
with an elliptical source is ten. 

3. SINGLE CAUSTIC CROSSING 

Due to the intrinsic difference in shape, the lightcurve of 
a caustic-crossing binary lens event involving an elliptical 
source differs from that of an event associated with a circu- 
lar source. However, detecting the deviation induced by the 
source flattening is difficult. This is in part because the flat- 
tening of the majority of stars is small. However, we find 
that even in the case of severely flattened stars the deviation 
caused by the source flattening is substantially smaller than 
the deviation caused by the limb darkening. In this section, 
we demonstrate this fact. 

To see the pattern of deviations induced by the source flat- 
tening and compare the magnitude of the deviation to that 
caused by limb darkening, we simulate events associated with 
source stars having various shapes, sizes, and limb-darkening 
coefficients. For the cases of events with circular sources, we 
consider limb darkening. We test four cases of limb darkening 
coefficients of = 0.1, 0.3, 0.5, and 0.7, which cover most 
values ranging from main-sequence to giant stars. The tested 
value of the normalized source radius is p^, = 0.02, which 
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roughly corresponds a giant source star of a typical Galac- 
tic bulge event. For events with elliptical sources, however, 
we do not consider limb-darkening effect, so as to isolate the 
magnitude of deviation caused only by source flattening. We 
test three different elliptical cases where the individual source 
stars have a common ellipticity but different orientation an- 
gles. For the elUpticity, we adopt a relatively large value of 
e = 0.5 {b/a = 0.866). The tested values of the orientation 
angle are = 0°, 30°, and 90°. We set the normalized semi- 
major axes. Pa = a/ 9^, of the individual elliptical source stars 
with different orientation angles so that the times of the caus- 
tic entrance and exit (and thus the duration of caustic crossing) 
match those of the circular source case. As a result, the indi- 
vidual elliptical source stars have different semi-major axes 
with ratios of paiijj = 0°) : pa(ip = 30°) : pai^p = 60°) = 1 .000 : 
1.035 : 1.110 (see the top panel of Fig.|5}- The duration of 
caustic crossing is proportional to the projected width of the 
source star in the direction perpendicular to the caustic, p± 
(projected width, see Fig.Q. The projected width is related 
to the semi-major axis and orientation angle by 



Ap_L = Pa 



(l-e^)2 + cot^V 



\-e^ + cor -0 



1/2 



(5) 



All lightcurves are produced by the same binary lens sys- 
tem with q = 0.8 and 5 = 1.2 and the parameters related to 
the source trajectory are mq = 0.27 and a = 0°. The angle of 
the source trajectory with respect to the perpendicular to the 
caustic, K (caustic incidence angle, see Fig.^, is very close 
to zero. 

In Figure 12] we present the results of the simulation where 
the middle panel shows the magnifications of the individ- 
ual tested events. A, and the bottom panel shows the devia- 
tions from a uniform circular case Aq. The top panel shows 
the shapes and orientations of the tested source star. The 
lightcurves of the limb-darkened circular source cases show a 
characteristic pattern in which the deviation is negative when 
the caustic is located on the fainter edge of the source while 
the deviation is positive when the caustic is positioned on the 
brighter central part. The elliptical cases show more com- 
plex pattern depending on the the size and orientation of the 
source star However, regardless of the variation, the devia- 
tion induced by the source flattening is substantially smaller 
than the deviation caused by the limb-darkening effect. 

4. MULTIPLE CAUSTIC CROSSINGS 

4.1. Basic Scheme of the Method 

Although difficult with only a single caustic-crossing 
lightcurve, it is possible to notice the elliptical nature of the 
source star if both the entrance and exit of the caustic are 
well resolved. This is possible because the products of the 
caustic-crossing timescale and the cosine of the caustic inci- 
dence angle, i.e. Af cosk, of the individual caustic crossings 
are different each other for an elliptical source. For a circular 
source, on the other hand, the products are always the same re- 
gardless of the source orientation. The product corresponds to 
the caustic-crossing timescale when the incidence angle of the 
source trajectory is k = 0° (hereafter we refer to the product as 
the 'normal incidence caustic-crossing timescale' and denote 
as /S.t±). The normal incidence caustic-crossing timescale is 
related to the source parameters by 



Af I = Af cos K = 2 



{l-e^Y + coV-ij; 
l-e^ + cot^ ip 
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Fig. 2. — Lightcurves of caustic-crossing binary lens events associated 
with source stars having various shapes, sizes, and limb-darkening coeffi- 
cients (middle pane!) and the deviations froin a uniforin circular case (bottom 
panel). Top panel shows the shapes and orientations of tested source stars. 



We note that the caustic-crossing timescale Af is defined as 
the duration between the source's enter and exit of the caus- 
tic. For each caustic crossing, the caustic-crossing timescale 
is measured from the lightcurve during the caustic crossing. 
The incidence angles k of the individual caustic crossings and 
the Einstein timescale Te are determined from the model fit to 
the overall shape of the lightcurve. 

Although the elliptical nature of the source star can be no- 
ticed from two caustic crossings, unique determination of the 
source ellipticity requires an additional crossing. The nor- 
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Fig. 3. — Relation between the source orientation angles ipi and ^2 with 
respect to the perpendiculars to first and second caustics, fti and K2 are the 
incidence angles of the source star to the first and second caustics, respec- 
tively. 



mal incidence caustic-crossing timescale is a function of three 
unknowns of pa, e, and tp, i.e. At± = At±(pa,e,tp). If one 
resolved the Ughtcurve profiles of two caustic crossings, the 
number of measured quantities increases to two, i.e. the nor- 
mal incidence caustic-crossing timescales at the first and sec- 
ond caustic crossings, Af 1 and At±,2, but the number of 
unknowns remains the same. This is because the source ori- 
entation angles of the individual caustic crossings are related 
each other. This is demonstrated in Figure |3l In the pre- 
sented case, the orientation angle at the second caustic cross- 
ing ip2 is related to the orientation angle at the first crossing 
tpi by ■02 = V"! + + ^^2, where ki and K2 are the caustic in- 
cidence angles at the first and second caustic crossings, re- 
spectively. However, with two caustic crossings, the number 
of unknowns still exceeds that of that of measured quantities. 
Therefore, if one has an additional lightcurve of caustic cross- 
ing, the numbers of unknowns and measured quantities be- 
come the same, and thus one can uniquely determine all three 
unknowns of pa, e, and ■0i (and subsequently ip2 and t/j^). 

There are two channels for multiple caustic crossings. One 
is the case in which the source trajectory passes two sepa- 
rate sets of caustics. The other channel is the case where the 
source trajectory is aligned closely along a fold of the caus- 
tic as shown in an example in Figure |3 Multiple crossings 
through the second channel can occur because fold caustics 
have concave curvature. If the individual caustic crossings 
of these events are resolved, the events will provide precious 
chances of measuring the source ellipticities. 

4.2. Feasibility of the Method 

The basic requirement of the proposed method is resolving 
multiple caustic crossings. Currently, microlensing follow-up 
observation is possible only after an alert is issued. Under this 
strategy, the first caustic crossing is usually missed due to the 
unpredictable nature and short duration of the caustic cross- 
ing. The second caustic crossing may be resolved because the 
crossing can be roughly predicted from the characteristic "U" 
shape of the lightcurve when the source is inside the caustic 
(jjaroszvhski & Mao 2001). In addition, there is a guarantee 
that the second crossing should occur and thus it is worth de- 
voting observational resources for the resolution of the second 
caustic crossing. However, there is no characteristic pattern 
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Fig. 4. — An example event with multiple caustic crossings. The small 
ellipses on the source trajectory (straight line with an arrow) represent the 
source star at the times of the caustic crossings, which are numbered ac- 
cording to the order of time. The size and flattening of the source star are 
exaggerated to better show the orientations of the major axis. The angle 9 is 
the orientation angle of the source star major axis with respect to the binary 
axis. The small bottom panel shows the lens system geometry, where the po- 
sitions of the lens components (mj and 1112), caustics, and source trajectory 
are marked. 



between the pairs of the caustic crossings, and thus it would 
be difficult to predict when the additional crossings will occur 
Furthermore, there is no guarantee that these additional cross- 
ings will occur. Therefore, considering that the main object 
of the follow-up observation is detecting deviations induced 
by planetary companions, it would be difficult to devote the 
limited resources waiting for unpredictable additional caustic 
crossings. 

However, the limitation can be overcome with the advent of 
future lensing experiments that can survey wide fields contin- 
uously at high cadence by using very large format cameras. 
Two such surveys were proposed. The Galactic Exoplanet 
Survey Telescope (GEST), whose concept was succeeded by 
the Microlensing Planet Finder {MPF), is a space missio n to 
be eq uipped with a 1 .Om - 1 .5m telescope ( Bennett & Rhid 
l2002h . The 'Earth Hunter' project plans to achieve high mon- 
itoring frequency by using three 2m-class wide-field ground- 
based telescopes scattered over the southern hemisphere (A. 
Gould, private communication). These next-generation sur- 
veys will dispense with the alert/follow-up mode and instead 
simultaneously obtain densely sampled lightcurves of all mi- 
crolensing events in the field. These experiments are expected 
to detect of order 5000 events per year Therefore, the set of 
multiple caustic-crossing events for which an ellipticity mea- 
surement is possible would not be vanishingly small. 

To investigate the feasibility of implementing the proposed 
method in future lensing surveys, we estimate the magni- 
tude of the difference in the normal incidence caustic-crossing 
timescales between caustic crossings, At±j- At±j, expected 
from typical Galactic binary lens events. As an example, we 
choose an event produced by a binary lens system with a mass 
ratio and separation of ^ = 0.1 and s = 1.8, respectively. Fig- 
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TABLE 1 

Normal Incidence Caustic-crossing Timescales 



e 


crossing 








Afx (hr) 




max|Arx , - Ar^ / I (hr) 


n 


number 


(°) 


(°) 


e = 0.2 


e = 0.3 


e = OA 


e = 0.2 e = 0.3 e = 0.4 



30 


1 


9.718 


12.718 


1.9182 


1.9162 


1.9137 


0.0359 


0.0820 


0.1488 




2 


77.689 


80.689 


1.8823 


1.8342 


1.7649 










3 


80.424 


77.424 


1.8831 


1.8363 


1.7691 










4 


53.208 


50.208 


1.8976 


1.8703 


1.8339 








60 


1 


9.718 


42.718 


1.9026 


1.8819 


1.8548 


0.0293 


0.0667 


0.1204 




2 


77.689 


69.311 


1.8863 


1.8438 


1.7839 










3 


80.424 


47.424 


1.8995 


1.8746 


1.8418 










4 


53.208 


20.208 


1.9156 


1.9105 


1.9042 








90 


1 


9.718 


72.718 


1.8848 


1.8403 


1.7769 


0.0341 


0.0775 


0.1393 




2 


77.689 


39.311 


1.9049 


1.8870 


1.8638 










3 


80.424 


17.424 


1.9167 


1.9129 


1.9082 










4 


53.208 


9.792 


1.9189 


1.9177 


1.9162 









urel^shows the caustics and source trajectory of the example 
event. For the source star, we assume a normalized semi- 
major axis of Pa = 0.002, which corresponds to that of a FO 
main-sequence source star of a typical Galactic bulge event 
produced by a low-mass stellar lens with distances to the lens 
and source of Di ^ 6 kpc and Ds ^ E kpc. We test three dif- 
ferent ellipticities of e = 0.2, 0.3, and 0.4, which correspond to 
the axis ratios of b/a = 0.979, 0.954, and 0.917, respectively. 
We test three cases of source orientations where the angles 
between the source major axis and the binary axis are 9 = 30°, 
60°, and 90° (see the definition of 9 in Fig.©. The small el- 
lipses on the trajectory in Figure |3 represent the source star 
at the times of caustic crossings. We designate the individ- 
ual crossings by numbers according to the order of time. We 
note that the source star in the figure is bigger and flatter than 
the actual tested stars to better show the orientations of their 
major axes. 

The results are summarized in Table [l] where we present 
the crossing number, caustic incidence angle k, source ori- 
entation angle ip, and the normal incidence caustic-crossing 
timescales Af^ of the individual caustic crossings. Also 
presented are the maximum values of the differences be- 
tween the normal incidence caustic-crossing timescales, i.e. 
max|Af^., - Af^j|. The range of the difference is order of 
minutes depending on the source orientation and ellipticity. 

Then, what would be the precision of At measurement? In 
an ideal case where there are equally spaced measurements 
over the caustic crossing with a fractional photometric preci- 
sion (7, the times of the caustic entrance and exit can be known 
with a precision of ^ a At /\/N. In reality, however, the actual 
error, 5{At) would be larger than this due to various factors. 
Estimating error in At measurement considering all these fac- 
tors is difficult. We, therefore, estimate the ratio of the actual 
error to the value of the ideal case, i.e. / = 6{At) / (a At / \/N), 
based on the error estimate in actually observed caustic cross- 
ings. For this, we choose two events wit h well resolved caus- 
tic crossings of MACHO 98-SMC-l (Albrow et al.' 1999b') 
and OGLE-1999-BUL-23 (Albrowetal. 2001). Assuming 
a = 2%, we find that / ~ 2.7 and 1 .9 for the individual events. 
By taking / = 2.0 as a representative value, we then estimate 



the precision of At measurement in the next-generation ex- 
periments as 

5iAt)^f-^. il) 



By assuming a typical caustic crossing timescale of Af ^ 2 
hrs, the photometric precision of a ^ 0.5%, and the monitor- 
ing frequency of 6 times/hr (and thus 12), we find that 
6(At) ~ 0.3 min. Considering that the typical value of the 
difference \ At±^i- At± j\ is order of minutes, we predict that 
ellipticity measurements will be possible for a significant frac- 
tion of multiple caustic-crossing events involving source stars 
having non-negligible ellipticities. 

5. CONCLUSION 

We have proposed a method that can determine the ellip- 
ticities of source stars of microlensing events produced by bi- 
nary lenses. The method is based on the fact that the normal 
incidence caustic-crossing timescales of the individual caus- 
tic crossings are different for events involving an elliptical 
source, while the timescales are the same for events associated 
with a circular source. For the unique determination of the 
source ellipticity, resolutions of three and more caustic cross- 
ings are required. Although this requirement is difficult to 
achieve under the current observational setup, it will be pos- 
sible with the advent of future lensing experiments that will 
survey wide fields continuously at high cadence. For typical 
Galactic bulge events, the difference in Afj^ between caustic 
crossings is order of minutes. Considering the high monitor- 
ing frequency of the future lensing surveys and the improved 
photometry especially of the proposed space-based survey, we 
predict that ellipticity determinations by the proposed method 
will be possible for a non-negligible number of events. 
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